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We consider a renormalization group study of the problem of coupled stripes of holes in cuprates. 
We use a model of a mesh of horizontal and vertical stripes and study the problem of superconduc- 
tivity via the Josephson coupling. We discuss the evolution of the Fermi surface with doping and 
temperature, the existence of Luttinger and/or Fermi liquid behavior, the presence of pre-formed 
Cooper pairs and the symmetry and magnitude of the superconducting order parameter. 
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After more than ten or so years after their discovery, 
there is still a debate about the microscopic physics of 
high temperature superconductors. Due to strong elec- 
tronic correlations transition metal oxides show a rich 
variety of phenomena such as anomalous transport, in- 
commensurate magnetic fluctuations and phase separa- 
tion. Incommensurate magnetic effects have been ob- 
served in La2-xSr x Cu04 and have been shown to 
have their origin in the formation of striped phases in re- 
lated compounds such as Lai.6-xNdo.4Sr x Cu04 0] and 
antiferromagnetic insulators such as La2- :E Sr. E Ni04+, / || . 
Recently it has been reported that the same incom- 
mensurate effects are observed in YBa2Cu307_a and 
Bi2Sr2CaCu208-a; Furthermore, various different ex- 
periments point towards a low dimensional picture of 
these compounds and to a closeness to quantum criti- 
cal behavior which can be interpreted in terms of static 
stripe formation ||. In the stripe picture the low di- 
mensionality is due to charge segregation into domain 
walls ||. Incommensurability appears because of the 
magnetic domain formation. In the antiferromagnetic 
phase stripes seem to explain well the available data . 
In this paper we study the problem of superconductivity 
and Fermi surface evolution with doping and temperature 
in a model proposed recently in which the stripes form 
a mesh of horizontal and vertical Luttinger liquids that 
hybridize with each other |H|. A similar model was pro- 
posed sometime ago in order to study the crossover from 
repulsively interacting one-dimensional Luttinger liquids 
to two-dimensional Fermi liquid s pA and has similarities 
with models of c-axis tunneling [[KJ. 

Our analysis is based on renormalization group calcu- 
lations ||il|] and we show that the model not only provides 
a very simple picture of the Fermi surface in underdoped 
cuprates but also agrees with recent experiments, espe- 
cially angle resolved photo-emission (ARPES) [p^|-p^[. 
We also calculate the relevant energy scales for the model 
and propose new experiments that can test our ideas. 
Our model differs substantially from earlier work on the 
effects of stripes in ARPES fb|. 

The model consists of horizontal and vertical stripes 
in different Cu02 planes described as Luttinger liquids 
with a characteristic interaction parameter g where g > 1 



(g < 1) for attractive (repulsive) interactions []T6| and 
separated from each other by antiferromagnetic regions. 
The Hamiltonian of the problem is written as (we use 
units such that h — ks = 1) 



TIt 



where Ti xn describes a ID chain with lattice spacing 
a in the x direction, which intersects the y-axis at n y £ 
where n y is an integer and £ = Na is the inter-stripe dis- 
tance. A similar notation is used to define TC Vi n w ■ TLt de- 
scribes tunneling between vertical and horizontal chains 
with amplitude t. We assume that the coupling is weak 
and treat TIt as a perturbation on the Luttinger liquid. 
We can define an extended Brillouin zone, in terms of 
the lattice constant of each individual chain, a. Each set 
of chains gives rise to a band in the extended zone. In 
the absence of hybridization the electrons are localized 
on the stripes and the Fermi surface is shown in Fig{j](a) 
where the horizontal and vertical bands are bounded by 
the lines at ±/c_f where kp = im/2 is the Fermi momen- 
tum for a chain with linear density n. 

The scaling equations for the coupling constants are 
obtained by tracing out high energy degrees and rescaling 
towards low energies [jll). The scaling equation for the 
tunneling is given by 

dit = C4I 

where £4 = 1 — 1/g where 1/g = (g + 1/g) /A and t is the 
ratio between the hopping in TLt and the infrared cutoff, 
A (all = cZ| A I / A) ]l7j . The renormalized tunneling energy 
is t R = t(t/w c )( 1 - C4 )/ £4 (a/^) 1 / < = 4 where uj c a bandwidth 
cut-off U . In the absence of tunneling, the electrons have 
a well defined momentum in the direction parallel to the 
stripe where they are located. TLt induces the folding of 
the Brillouin zone shown in Fig.[j](a) into N x N smaller 
zones. However, these are higher order effects in terms of 
the tunneling and distortions of the mesh, induced either 
by static disorder or thermal fluctuations |l8| , will smear 
these additional "shadow" bands |l9| . Inter-chain tunnel- 
ing will also hybridize the states derived from horizontal 
and vertical bands. Observe that for a doping amount of 
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x holes there are 1 /N — x electrons participating on this 
Fermi surface. The other 1 — 1 /N electrons form the sur- 
rouding antiferromagnetic background. Since the antifer- 
romagnetic and stripe electrons are in equilibrium they 
must have the same chemical potential (so that the lower 
Hubbard band of the antiferromagnet is pinned at the 
stripe Fermi surface). Thus, superimposed to the Fermi 
surface of Fig.|l|(a), one must have a diamond shaped, and 
somewhat blurred, occupation in momentum space for 
the antiferromagnetic electrons, with area proportional 
to 1 - 1 /N. The final result is shown on Fig.|l|(b). The to- 
tal area seen by ARPES should be proportional to 1 — x, 
corresponding to a large Fermi surface. These results 
agree well with the available ARPES data [p|,|3). 
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FIG. 1. (a) Fermi surface for the stripes without hybridiza- 
tion; (b) Fermi surface with hybridization: the continuous line 
shows the region of Fermi liquid behavior (finite residue) and 
the doted line shows region of non-Fermi liquid behavior. 

In most points of the Brillouin zone, the horizontal 
and the vertical bands are not degenerate and there is 
a splitting Sk m ,k y between bands. At zero temperature 
<5fc x ,fc y cuts off the RG flow. If 5k x ,k y > tR the verti- 
cal and horizontal stripes decouple and coherence effects 
are not expected. As shown in ref. |1 the quasiparti- 
cle residue vanishes in the regions close to (±rr/a, 0) and 
(0, ±7r/a) and the system behaves like a Luttinger liq- 
uid. The only exception are the points close (±fcp,±fcir) 
where S±k F ,±k F = 0. When Sk x ,k y < £r the horizontal 
and vertical bands will become strongly hybridized, and 
a quasiparticle residue is finite [^o| leading to local Fermi 
liquid behavior. The crossover from Luttinger to Fermi 
liquid along the Fermi surface is defined by tR f=a <5fe x ,fc H as 
shown in Fig.[l](b) . It turns out however that the regions 
of non- Fermi liquid behavior are nested and therefore will 
be subjected to instabilities pi] ]. An interesting conse- 
quence of this scenario is that there is a crossover from 
Luttinger liquid to Fermi liquid behavior also as a func- 
tion of £ as £ decreases. There is a minimum value of £, 



say £*, for which S(£*) w tail*). Since £ is a monontonic 
decreasing function of x B there is a crossover doping x* 
such that for x < x* (x > x*) we have Luttinger (Fermi) 
liquid behavior. 

Interaction processes do not need to involve an energy 
cost, even if they conserve crystal momentum. We can 
use a standard diagrammatic analysis to study the in- 
teractions which can be generated by the combined ef- 
fect of the intra-chain couplings and the inter-chain hop- 
ping. The irreducible diagrams that can be built starting 
with two electrons in either the horizontal or the vertical 
chains, which interact between them, hop into the other 
set of chains, and interact again. Besides the usual pair- 
ing process within each chain which is shown in Figj^(a) 
there are two types of possible couplings which satisfy 
momentum conservation which are depicted in Fig.^(b) : 

i) The Cooper pair channel: Cooper pair in the horizontal 
chains is transferred to the vertical chains, or vice- versa; 

ii) The transfer of Cooper pairs with total momentum 
(±2kF, ±2/cf). There are also direct (Fig.||(c)) and ex- 
change (Fig.|2j(d)) diagrams which give a vanishing con- 
tribution ibecause of phase space restrictions jllj . These 
processes can all be written as A v h — 0\Oh, where O 
is an operator which creates, or destroys, a pair of elec- 
trons in one chain. As each separated chain is a Luttinger 
liquid, the O's acquire anomalous dimensions relative to 
the non-interacting fixed point. The successive applica- 
tion of these operators generate interactions which take a 
pair of electrons from one set of chains into the other set 
and back again. Thus, upon scaling towards low energies, 
interactions of the type A^h and A vv will also be gener- 
ated. In general, the operator A has to be considered a 
2x2 matrix. It can be shown ||| that the O associated 
with processes ii) is less relevant than the operator which 
describes processes of type i). Note that, as our scheme 
implies a summation to all orders of intra-chain effects, 
the propagators shown in Fig.|| acquire anomalous di- 
mensions ]lq| . The scaling equation for the matrix Aij 
has the Riccati form: 

8 t A = -cx{g, l)Pa x + 2c 2 (g)A - c 3 (g, l)A 2 (1) 

where a x is a Pauli matrix. c\(g,l) — cie 4 *- 9-1 -*' repre- 
sents the generation of the Cooper pair tunneling start- 
ing from the tunneling of individual electrons. In a con- 
ventional junction between superconductors, this is the 
Josephson coupling (ci > unless the tunneling takes 
place through a magnetic region |2l| ). The second term 
is the scaling dimension of the operators with ci > 
(ca < 0) for attractive (repulsive) intra-chain interac- 
tions. c^{g,l) = e 2 *- 9-1 -*' is associated with the two- 
dimensional pairing on the Fermi surface in Figj^(a) [ pd| . 
Observe that due to symmetry A = AI + Ba x where A 
and B are real numbers and / is the identity matrix. It 
is clear from equation (|l]) that it is useful to work with 
the eigenvalues of A: 

\± = A±B. (2) 
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After a redefinition of the variables (A± 
r = e 3 ^ 9_1 ^f) the RG equations read: 

a i A+ = -c 1 r 2 + 2c 2 ( 5 )A + 
d t \- =c lT 2 + 2c 2 (g)\- - A 2 
diT = c 4 (g)r 



,2(9-1)1 



A±, 



A 



(3) 



where 02(g) — g - 1 + c 2 and 04(g) — 3(g — 1) + 24(5). 
The RG flow has four fixed points: (0,0,0), (2c 2 ,0,0), 
(0,2c 2 ,0) and (2c 2 ,2c 2 ,0). The (0,0,0) is the non- 
interacting fixed point; (2c 2 , 2c 2 , 0) is the non-interacting 
fixed point of isolated chains but with superconducting 
fluctuations. The physics of the other two fixed points 
is more complex. To gain more insight we study the RG 
flow when r = 0. 



of repulsive interactions the mesh should be insulating 
because of the intrinsic disorder in the Cu0 2 planes [0 . 
A mesh like the one proposed here is very sensitive to dis- 
order and impurities, such as Zn. Note that Zn reduces 
drastically the superfluid density |2q]. The attraction 
within each stripe is provided by the surrounding anti- 
ferromagnet |j|,^6| . The model proposed here should also 
be useful in describing artificial Josephson arrays pTjj . 

Since the flow is divergent at some finite scale /* we 
estimate the gap from the scaling equations by studying 
the flow close to the singularity. At this point the tun- 
neling amplitude is constant and given by r» = t(/*). If 
T * > c 2 /ci the gap has a BCS form, 
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FIG. 2. (a) Cooper pair diagram within a horizontal stripe 
(H); (b) Cooper pair diagram for a pair in a horizontal stripe 
(H) tunneling (T) to a vertical stripe (V); (c) direct processes; 
(d) exchange processes. 

For attractive interactions the flow is shown in 
Fig.||(a). Observe that (2c 2 ,2c 2 ,0) is stable whereas the 
non-interacting fixed point is unstable; the other two 
fixed points are also unstable. (2c 2 , 0, 0) is unstable in the 
negative A + direction and from the diagram of Fig.^|(b) 
we see that the instabilities is of the superconducting 
type. One can understand the strong coupling fixed point 
in a very simple way: observe that A_ is constant and 
A+ << -1. From (0) we have A < B « -1. It implies 
a superconducting instability in which the sign of the or- 
der parameter is the same along horizontal and vertical 
stripes; this is equivalent to s-wave pairing. (0,2c 2 ,0) is 
unstable in the negative A_ direction while A+ is con- 
stant which implies A « —I and B » 1. The co- 
efficients have opposite sign which means that the sign 
of the order parameter in the vertical stripes is opposite 
to the sign in the horizontal stripes leading to what is 
called a it -junction fl23fl . Since the stripes are oriented 
along the crystallographic directions this is equivalent to 
d-wave pairing. For repulsive interactions the s-wave and 
d-wave pairing persist and the flow is shown on Fig|j|(b). 
If r 7^ the hopping is relevant (except for strongly at- 
tractive interactions). If c\ > the system flows towards 
s-wave and if c\ < the flow is reversed towards d-wave 
superconductivity. One also expects that in the presence 



where Ao is the coupling constant and u> p is the charac- 
teristic pairing energy. In this limit one expects usual 
BCS behavior with well defined quasiparticles in the su- 
perconducting state. Eq. (0) reduces to the usual BSC 



gap (Abcs 



-l/IAol 



) in the non-interacting limit 



(ci, 2 — > 0). In the opposite limit, 77 < c 2 /ci, Luttingcr 
liquid behavior is obtained with a non-BCS gap given by 



A - \fc 



C-2 



A + \f(- 



C-2 



which, when c\ — > 0, has the interesting form = 
u}p(l + c 2 /| Ao|) -1 / C2 which raises much faster with the 
coupling constant than the usual BCS result. 

Let us consider the tunneling of Cooper pairs between 
transverse stripes A and B through the antifcrromagnetic 
regions. Since the pair is in a singlet state the tunnel- 
ing process occurs in such a way that a given spin or- 
der the sign of the wavefunction is reversed ^3|. This 
is the so-called a 7r-junction and implies that c\ < in 
the RG flow. The argument is valid if the tunneling oc- 
curs through a re gion smaller than the magnetic corre- 
lation length, £jvf IPg| - A similar mechanism for d-wave 
pairing involving tunneling through a finite antifcrromag- 
netic region was found by Bonesteel (29) . Assuming that 
ci < the tunneling through the antiferromagnet only 
takes place if £m{T*) ~ a. Thus, there is a crossover 
temperature T* at above which c\ = [Q. T* can be 
roughly estimated assuming that the antiferromagnet is 
in the quantum critical regime p8[ | where T* w c/a (c is 
the spin wave velocity). 
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FIG. 3. Renormalization group flow in the t — plane: (a) 
attractive interactions; (b) repulsive interactions. 



ior) and at lower temperatures a spin gap appears before 
the system superconducts. 
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Note: Since completing this work, we became aware 
of numerical work on stripes in the t-J model ]33| which 
seems to be consistent with our results. 



In conclusion, besides the shape of the Fermi sur- 
face, the existence of "shadow" bands, the co-existence 
of Fermi liquid and non- Fermi liquid behavior, the pos- 
sibility of non-BCS and BCS behavior and s- and/or d- 
wave pairing we have three different temperature scales 
in this problem, namely, T* , S and A. Consider first 
the case such that T* > S > A. If T* > T > S the 
pairing takes place but there is no coherence between 
the stripes. There are pre-formed pairs and a gap to 
spin excitations due to the formation of singlets. This 
effect is observed experimentally in infrared absorption 
and nuclear magnetic resonance (NMR) |n|. If 6 > T 
tunneling takes place and a gap opens at the (±7t/a,0), 
(0,±7t/a) regions of the Fermi surface but no phase co- 
herence is observed. This case has similarities to a recent 
model proposed by Geshkcnbcin et al. for the physics 
of the Fermi surface with pre-formed pairs J3^] . In the 
other scenario, S > T* > A, the physics is quite dif- 
ferent. For S > T > T* a Luttinger liquid is formed 
close to (±7r/a, 0), (0, ±n/a) but no Cooper pairs are 
found. When T* > T the Luttinger liquid regions of the 
Fermi surface become gapped due to pair formation. At 
lower temperatures superconductivity is obtained. While 
the first case is closer to a conventional superconducting 
instability the last one shows direct transition from a 
Luttinger liquid to a superconductor fiofl . This picture 
agrees with the available data in the underdoped cuprates 
p2| , ^3| and provides a microscopic origin for the boson 
modes proposed by Shen and Schrieffer jl4j in terms of 
collective modes of the Luttinger liquid. Observe that in 
the underdoped compounds we predict weak dependence 
of Jcf on doping [ p8j . 

Finally, we can predict a new effect as we change the ra- 
tio r w T*/6 « (c/vp)(£/a). Close to the metal- insulator 
transition one expects r > > 1 because c can be large and 
Ija w l/x » 1. In this case the first scenario applies 
and as the system is cooled down a spin gap must be 
observed before a gap opens at (±n/a, 0), (0,±n/a) in 
the Fermi surface. At larger doping (but smaller than 
x*), that is, when r << 1, one expects that no gap will 
appear in the Fermi surface (because of Luttinger behav- 
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